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ABSTRACT 

Context. Galaxies are believed to be the main providers of Lyman continuum (LyC) photons during the early phases of the cosmic 
reionization. Little is known however, when it comes to escape fractions and the mechanisms behind the leakage. To learn more, 
one may look at local objects, but so far only one low-z galaxy has shown any signs of emitting LyC radiation. With data from the 
Far Ultraviolet Spectroscopic Explorer (FUSE), we previously found an absolute escape fraction of ionizing photons (f esc ) of 4-10% 
for the blue compact galaxy Haro 1 1 . However, using a revised version of the reduction pipeline on the same data set, Grimes and 
collaborators were unable to confirm this and derived an upper limit of f esc <2% . 

Aims. We attempt to determine whether Haro 1 1 is emitting ionizing radiation to a significant level or not. We also investigate the 
performance of the reduction pipeline for faint targets such as Haro 11, and introduce a new approach to the background subtraction. 
Methods. The final version of the reduction pipeline, CalFUSE v3.2, was applied to the same Haro 11 data set as the two previous 
authors used. At these faint flux levels, both FUSE and CalFUSE are pushed to their limits, and a detailed analysis was undertaken 
to monitor the performance of the pipeline. We show that non-simultaneous background estimates are insuffient when working with 
data of low signal-to-noise ratio (S/N), and a new background model was developed based on a direct fit to the detector response. 
Results. We find that one has to be very careful when using CalFUSE v3.2 on low S/N data, and especially when dealing with sources 
where signal might originate from off-center regions. Applying the new background fit, a significant signal is detected in the LyC in 
both detector segments covering these wavelengths. Thus, the leakage is confirmed with a flux density of /9oo=4.0 x 10~ 15 erg s -1 
cm -2 A" 1 (S/N=4.6), measured on the airglow free regions in the LyC for the night-only data. This corresponds to an absolute escape 
fraction of ionizing photons from Haro 1 1 of f esc =33±0.7 %. We confirm these results by investigating the two-dimensional data, the 
count rates, and the residual flux in Cn /11036A. 

Key words. Galaxies: intergalactic medium - starburst - fundamental parameters - evolution, Cosmology: diffuse radiation, 
Ultraviolet: galaxies 



1. Introduction 

One of the main events in the history of the universe was the 
phase transition when the neutral gas again became ionized. 
Finding out when and how this cosmic reionization occured 
remains a challenge for observers in a broad range of fields. 
Observational constraints based on the seven-year WMAP data 
indicate that the reionization started at a redshift z»T0.5 ( |Larson| 
|et al.||20T0 i. Over time, the intergalactic medium became ion- 
ized over larger and larger regions, but the redshift at which this 
process was completed is uncertain. While some studies of the 
Gunn-Peterson effect, the evolution of the Lya escape fraction, 
and the Lya luminosity function indicate that reionization ended 
at z^6 (e.g |Malhotra_& Rhoads|[2004l |Fan et al.|[2006l |Totani| 



et al.|200"6 Ono et al.|2010|l, other investigations give less sup 



port for this (e.g. Tilvi et al.|2010| |Mesinger|2010| Ouchi et al. 
|2010| l. A highly inhomogeneous intergalactic medium compli- 
cates the interpretation of the scanty observational data but if we 
assume that the reionization epoch took place over the redshift 
interval z=6-ll, we can try to identify the sources generating 
the cosmological background of ionizing photons. Probably the 
most significant sources varied over time. At late times, z<3, it 
has been found that the bulk of the ionizing photons originates 



from quasars ( |Fauch er-Gigu ere et al.| |2008) and AGNs (Cowie 
|et al.|2009] >, while at higher redshifts quasars are too few in num- 
ber to be the dominant source ( |Jiang et aT 2008 Loeb|2 009 ). The 
onset of the reionization has instead been connected to the ap- 
pearance of the first luminous structures in the universe, the first 
stars and galaxies piaehnelt et al.|200T 



It still remains unclear whether there were enough galaxies 
around at the epoch of reionization to account for the required 
number of ionizing photons, but lately several high-z surveys 
provide support for an early galaxy build-up. Using the HST 



WFC3/IR camera, |Bouwens et al. ( 2009} derived the galaxy lu- 
minosity function for z « 10 galaxies, and found an increasing 
likelihood of a galaxy-driven reionization. The evolution of the 
luminosity function together with models show that the onset 
of reionization must have been initiated by low mass galaxies 
(e.g. |Trenti et al.|2010[|Srbinovsky & Wyithe|2010| , with prop- 
erties similar to the observed low luminosity population at z~7 
(Bou wens et al.|2010)>. Cosmological simulations i ncluding stel- 
lar feedback (e.g. |Razoumov &~S ommer-Larsen 2010, Yajima 
|et al.|20iTj| give further support to the importance of LyC leak- 
age from dwarf galaxies (halo masses M/, « 10 9 Mq). The same 
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models show a decreasing contribution from galaxies of increas- 
ingly higher mass (but see Gnedin et al.|2008 1 



Direct measurements of escaping ionizing photons from the 
galaxies that reionized the universe is unfortunately impeded by 
the still highly opaque neutral intergalactic medium (IGM) at 
these redshifts. To measure escape fractions and investigate the 
physical properties that allow leakage, we need to study more 
nearby objects. This can and has been done with groundbased 
telescopes at z~3 (e.g.|Steidel et al |2001| |Shapley et al 12006} 
|Iwata et al.|p009) |Vanzella et al.||2010| |Ne"stor et al.||2011), 
or locally by using space telescopes suc h as HUT (|Leitherer| 
et al.|1995), FUSE (jPehirveng et al.|2001[|Bergvall et al.|2006| 



Grimes et al.|20 07 2009 ), GALEX (|Cowie et al.|2 009 ), and HST 



(Malkan et al. 2003, Siana et al. 2007 2010 Bridge et al. 2010). 
Of the four only HST and GALEX are still in operation, and un- 
til now only HUT and FUSE have been providing wavelength 
coverage down to 900 A, such that local galaxies can be stud- 
ied. With the successful installation of the new HST instrument 
COS, this wavelength region is yet again opened up. Despite all 
of these efforts, only a few investigations have resulted in LyC 
detections for single sources. There have been ~50 single de- 
tections so far, all except one at z~3, although foreground con- 
tamination is expected to account for some of those. The one 
exception is the disputed case of Haro 1 1 . 

Over the years, FUSE has provided a wealth of information 
about the ultraviolet (UV) spectrum of stars and galaxies and 
made an outstanding contribution to UV astronomy. It was how- 
ever not optimized for observing faint objects, and its sensitiv- 
ity is at its lowest at the wavelengths that are ideal for measur- 
ing the Lyman continuum from local galaxies. Therefore, even 
though the FUSE archive still provides a unique legacy to aid 
the search for local LyC leakers, the data reduction has turned 
out to be non-trivial and few results have been published on the 
subject. A first investigation of the local galaxy Mrk 54 resulted 
only in an upper limit to the absolute LyC escape fraction f esc < 
6% ( |Deharveng et al.||200"T) l. In |B erg vail eTal] ( j20"06| >, the first 
actual detection of LyC es cape w as reported for Haro 11 to be 
/ esc =4-10 %. However, in Grimes et al.| (2007), the same data 



set was reduced with a later version of the FUSE pipeline, and 
here they arrived at contradicting results with f esc < 2 %. The 
LyC escape was studied in eight more galaxies from the FUSE 
archive in Grimes et al.|(2009 ). None were found to have any sig- 



nificant excess in the LyC, and several display a negative flux in 
the relevant wavelength region as a result of the over-subtraction 
of the reduction pipeline background. 

In this article, we present a detailed analysis of the same 
Haro 1 1 data set as used by Bergvall and Grimes, to determine 
whether Haro 1 1 is a LyC leaking galaxy. In addition, we ana- 
lyze the final version of the FUSE reduction pipeline, CalFUSE 
v3.2, itself in terms of its treatment of faint targets at the short- 
est wavelengths. We introduce our own background fitting pro- 
cedure and an alternative way of extracting low signal-to-noise 
ratio (hereafter S/N) FUSE spectra. 



2. Haro 11 

Haro 11, also known as ESO 350-IG38, is a local blue com- 
pact galaxy with M B =-20.3. It seems to share many proper- 
ties with star-forming galaxies at z~3 (Lyman break galaxies, 
LBGs), although being at the bright end of the far-UV luminos- 
ity distribution (e.g. Overzier et al. 2008 ). Haro 1 1 has a complex 
morphology with three main star-forming knots and se ems to 
be dynamically unrelaxed, remniscent of a merger event ( Ostlin 




Fig. 1. Left: Composite HST image of Haro 11. Image scale is 
15x15". East is left, north is up. The Lya emission is shown 
in blue, the ionized gas via Ha in red and the UV-continuum 
in green. The three main star-forming knots A, B and C are la- 
belled. Right: The Lya vs. the Ha surface brightness for the 
binned elements in the image. The hatched lines show Lya=Ha 
and the case B Lya/Ha=8.7. At low surface brightness (1) the 
theoretical ratio is exceeded due to resonant scattering mak- 
ing the Lya emission off-set from the ionized regions, while 
at higher surface brightness (2) the top edge of the distribu- 
tion agrees with the recombination value, indicating that the Lya 
photons from these regions are scattered just a few times before 
escaping. Regions with heavy extinction can also be seen at high 
Ha surface brigh tnesses (3). The fig ure is an altered composite 
of Fig. 3 and 4 in |0sflin et al. (2009 1 (see this article for further 
information). 

(A color version of this figure is available in the online journal.) 



et al. 



200l|. The knots each contain numerous super star clus- 



ters, for which Adamo et al. ( 2010j ) derived young ages, between 
1 and 40 Myr and masses between 10 3 and 10 7 M0. Haro 11 is 
also found to have a low metallicity, with an oxygen abundance 
12+log(0/H)=7.9 (Ber gvall & 0stlin|2002t , and a surprisingly 
low H I content. The estimated upper hmit to its Hi mass was 
found to be M(H I) < 10 s M in |Bergvall et al.| ( |2006| l, who also 
noted that Haro 1 1 seems to contain very extended photodisso- 
ciation regions (with H2 molecules) of a total mass exceeding 
the upper H I mass limit. The same authors also found a similar 
amount of ionized gas in Haro 1 1 . Thus the dominating hydro- 
gen phases could be H + and H2, providing more favourable con- 
ditions for escape than in normal galaxies where H I serves as an 



effective shield for the LyC radiation. |Bergvallet al. (20061 also 
investigated the ionizing stellar population and found signatures 
of Wolf Rayet stars in the optical as well as P Cygni profiles 
in the far-UV, revealing the presence of O supergiants and an 
ongoing burst or a burst that has only very recently ended, in 



agreement with the ages derived in Adamo et al. (20101. 



Haro 11 is also known to be a Lya emitter (e.g. ( |Kunth 
|et al.| 1998 1, and the location where the Lya photons escape the 
galaxy has been established using high resol ution imaging with 
the ACS/SBC onb oard the HST spacecraft ( |0stlin et al.| |2009 



|Hayes et al.|2007 1. These studies found that the main part of the 



Lya photons are emitted in a diffuse halo after multiple scatter- 
ings in the interstellar medium, that is clearly offset from the ion- 
ized gas surrounding the star-forming regions (see Fig. [T] where 
the Lya is shown in blue, the ionized gas via Ha in red, and 
the UV-continuum in green). Of particular interest in the light of 
possible LyC leakage is one of the star-forming regions, knot C. 
The right panel in Fig. T] shows the Lya versus the Ha surface 
brightness for binned elements. For higher surface brightnesses 
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in knot C, it was found that the Lya /Ha ratio agrees with the the- 



oretical case B recombination value, Lya/Ha=8.7 (Osterbrock 



& Ferland|2006| l. This indicates that the Lya photons from this 
region are scattered just a few times before escaping. Such direct 
escape would be possible from a low-density, dust-poor nebula 
bordering on the optically thin case A, an environment that could 
also allow ionizing photons to escape. If the Lya photons es- 
caped by means of an outflowing neutral medium on the other 
hand, as observed in Kunth et al. (1998), such an escape mech- 
anism would not aid the LyC photon escape unless the wind is 
strong enough to clear channels. 

Knot C is certainly interesting in the light of possible leaking 
LyC radiation, but might not be the most likely option since the 
amount of neutral gas and dust should be largest in the central 
parts. Iwata et al. (2009 ) and Nestor et al. ( 201 1\ demonstrated 
that LyC leakage can be found spatially off-set from the main 
body of the non-ionizing UV continuum and in the central parts 
of the galaxies (although some of these detections statistically 
are likely to be lower redshift contaminants). To estimate the 
distance from the center at which we would expect most of the 
ionizing radiation to slip through a picket fence structure in Haro 
11, we adopt a simple model with two components - the ioniz- 
ing stars and the neutral gas. We would ideally wish to use a 
map of the UV flux, combined with a high resolution H i map. 
Since we do not have access to the se data we will instead u se 
the B band flux distribution T. B fi( r ) ( |Bergvall & 6stlin|2002 ) to 
model the ionizing flux, assuming a thin disk. Since H i has not 
been detected in Haro 11, we assume that the hydrogen gas is 
distributed in a homogeneous exponential disk with scalelength 
h.Hi- In a picket-fence model, the LyC photons are absorbed by 
spherical clouds with radius R and cross-section <x = jiR 2 . Thus 
the filling factor is a = n ■ ^f-, where n is the number of clouds 
per volume unit. The observed LyC flux is then 



fi 



L\C,obs 



= 2n f 
Jo 
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o(r) ■ re 
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where Sz, v c,o(r) is the surface brightness in the LyC (here as- 
sumed to be oc £ B X an d r is the distance from the centre. The 
optical depth, r(r) is 
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where z is the vertical height, and h HI ,- and h HIz are the radial 
and vertical scalelengths of the gaseous disk. 

We do not know the scalelength of the Hi disk but a rea- 
sonable guess is that it is roughly twice the optical scale length 
(Ivan Zee et al.pOOT] |Broeils & van_Woerden|| 19941 IBroeils &| 
|Rhee|1997| ), i.e. about 6 kpc ( jBergvall & Ostlin|2002| ). We reg- 
ulate cr until we obtain a leakage corresponding to 3%. We then 
find that the galaxy is leaking almost equal amounts at radii in 
the range 1-14 kpc and that ~50% of the radiation is leaking 
through outside ~1 kpc (or 17") from the center. This simple 
model indicates that Haro 11 is an extended source in the LyC 
as seen by FUSE, and that we cannot rely only on the height of 
the spectrum longwards of the LyC to make a qualified guess as 
to whether the source is extended. 

In Fig. [2] a Gunn r-band image of Haro 1 1 is displayed (thus 
including nebular Ha emission). A slit spectrum is shown to the 
right, taken at the position marked by the hatched lines in the im- 
age. The image was obtained with the Gunn r filter at the 1.5m 




Fig. 2. Haro 1 1 in a deep Gunn r image. East is left and north is 
up. The contour plots show the far-UV continuum at Lya, and 
help indicate the size of Fig.[T] The slit spectrum to the right was 
obtained at the position marked by the hatched lines in the im- 
age. A star-forming region is found located 35" from the center 
north of Haro 11, illustrating that these regions exists out to, and 
beyond, the 30"x30" field of view covered by the FUSE LWRS 
aperture. 



Danish telescope at La Silla (90 minutes, 1988) and the slit spec- 
trum with grism B300 using EFOSC1 at ESO 3.6m, La Silla 
(700s, 1984). In the central parts, contour plots of the far-UV 
continuum at L ya are shown, wh i ch have a size corre sponding to 
Fig.p](credit: |Hayes et al.| ( |2009| ; |6stlin et al.| ( |2009] )). The opti- 
cal size of Haro 1 1 is obviously much larger than the HST/ACS 
image shows. The slit spectrum reveals a highly blueshifted star- 
forming region at a distance of 9 kpc (35") . The stellar contin- 
uum in the optical is not visible, and the region is identified via 
its nebular emission. Although this object falls outside the FUSE 
aperture, it shows that star-forming regions do exist out to large 
distances where it might be easier for LyC photons to escape. 

The redshift of Haro 1 1 is z=0.021, which places the Lyman 
limit (Ly-limit) at 930.6 A. 



3. Observations, reductions, and data analysis 

The Far Ultraviolet Spectroscopic Explorer (FUSE) is a space 
telescope that operated between the years from 1999 to 2007. It 
was equipped with a far-UV spectrometer that covers the wave- 
length range 905-1187 A with a resolving power of ~ 20000. 
The four co-aligned telescopes each has a spectrograph built ac- 
cording to a Rowland circle design, where the spectra are im- 
aged on two microchannel plate detectors. To accomplish good 
throughput over the full wavelength range, the FUSE instrument 
was divided into complementary channels based on silicon car- 
bide (SiC) optics for wavelengths shorter than ~1000 A and 
lithium fluoride (LiF) over-coated aluminium for wavelengths 
in the range ~1000-1200 A. Each detector is divided into two 
segments that each image both a SiC and a LiF spectrum, so 
there are in total eight spectra. Two of these spectra cover the 
wavelength region where we wish to measure any possible LyC 
signal, SiC IB (905-992 A) and SiC 2A (916-1005 A). For de- 



4 



1 Elisabet Leitet et al.: Analyzing low signal-to-noise FUSE spectra: 



tails about the on-orbit performance of the FUSE instrument, we 
refer to |Sahnow etaTJ l |200Q} . 

Detailed information about the satellite and its instrument 
can be found in the FUSE Archival Instrument Handbook] and 
The FUSE Archival Data Handbook These two handbooks will 
be referred to as IH and DH in the text. 

The observations of Haro 1 1 were performed on October 12, 
2001, using the FUSE low resolution mode, LWRS, which has 
an aperture of 30"x30". The total exposure time was 16.5 ks, of 
which 12.2 ks was taken during orbital night. The OBSID of the 
data set is B 1090101. 



3.1. Reductions and data analysis 

3.1.1. The reduction pipeline, CalFUSE 

We now present a short introduction to the data reduction 
pipeline of FUSE, CalFUSE. CalFUSE has been under constant 
development since the start of the mission in 1999. Although 
the two previous articles on Haro 1 1 used earlier versions of the 
pipeline, here we only used the final version, CalFUSE v3.2, 
which is described in detail in Dix on et aT7| ( |2007| >. 

Running CalFUSE v3.2, the data are corrected for various 
motions and distortions introduced by the satellite and instru- 
ment. The arrival time and position of the photon events are 
preserved for time-tag observations (TTAG), and stored in two- 
dimensional (2D) intermediate data files (IDF). The photon list 
is kept until the last possible moment when the one-dimensional 
(ID) spectra are extracted, which makes it possible for the user 
to easily screen out data, e.g. for observations too close to the 
earth limb, during count rate bursts, during orbital day, or for 
photons corresponding to wavelength regions affected by air- 
glow. All photons are flagged by the pipeline, and the user can 
choose to filter out the unwanted ones. 

The background model of the pipeline is of special inter- 
est in this article. The model consists of three different parts, 
one spatially uniform originating from intrinsic 4() K decay in 
the microchannel plates together with cosmic rays, and two spa- 
tially varying day- and night components from scattered light. A 
detailed description of the internal detector background can be 
found in IH, Section 4.4.3.1 and of the stray and scattered light 
backgrounds in DH, Section 7.2.1. The pipeline first scales the 
intrinsic background and then fits fixed background files to the 
varying day- and night components. These template background 
files were produced by taking long exposures of a blank field 
of the sky repeatedly over the mission, and the two closest in 
time to the observation are selected and interpolated. The shape 
of the background is therefore pre-determined from those tem- 
plate files. This procedure works satisfactorily for bright targets, 
but for low S/N spectra the pipeline has been known to overes- 
timate the background at shorter wavelengths producing an un- 
physical negative flux (seen in e.g. |Grimes et al. 2007] |Grimes| 
et al.|2009| l. Particularly problematic are the SiC IB and SiC 1A 
LWRS spectra, because they fall close to the edge of the detec- 
tors where the background has a steep slope and is difficult to 
model properly. A 10 % uncertainty estimate in the background 
was made for earlier versions of the pipeline and propagated 
through to the final spectrum, but in later versions this is only 
treated as a systematic error and as such not included in the final 
errors ( |Dixon et al.|2007) >. 



In several previous works where FUSE data of low S/N were 
treated, non-standard procedures were applied in the data reduc- 
tion. While excluding orbital day data to minimize the scattered 
light component seems to be common in most of these investi- 
gations, they differ in other aspects. In Fechner et al. (2006]l, the 
pulse height range was modified to reduce the background, while 
in |Dixon et al.| ( [2006] l they kept the standard pulse height limits 
with the argument of preserving photometric accuracy. In Zheng 
|et al.| ( [2004] >, they extend the analysis and reduce the (night only) 
data with both the standard pipeline procedure and a manual ex- 
traction. In the manual method, two background regions were 
defined, one on each side of the spectrum, and the interpolation 
of the background between those made under the assumption 
of linear variation. For the spectra close to the detector edges, 
two regions were defined on one side of the spectrum and the 
background extrapolated from these. They found that the man- 
ual approach was more successful, as judged by the level of the 
residual flux in saturated absorption lines. 



3.1 .2. [Berg vail et al. (2006) and Grimes et al. (2007) 



1 http://archive.stsci.edu/fuse/ih.html 

2 http://archive.stsci.edu/fuse/dh.html 



There have been two previous papers published on the leakage 
of ionizing photons from Haro 1 1 based on the FUSE data. Here, 
a brief summary of the two is presented. 

In |Bergvall eTaLl ( [2006l ), the pipeline version CalFUSE v3.0 
was used with default parameters. Haro 11 was treated as an ex- 
tended source, and the data included both the orbital day and 
night observations. Ber gvall et al.| ( [2"0"0"6| ) detected a visible LyC 
signal in both the SiC IB and SiC 2A IDFs. The extraction of 
the ID spectrum by the pipeline showed a difference between 
the SiC IB and the SiC 2A flux density, and the authors cor- 
rected this by adapting the background level of SiC IB to the SiC 
2A. The LyC signal at rest wavelength 900 A was determined 
to be /0oo=Ll±O.lxlO erg crrT 2 s A -1 after correcting for 
Milky Way extinction. This gave an absolute escape fraction of 
/esc=4— 10%, depending on the choice of parameters to describe 
the spectral evolutionary distribution. The authors pointed out, 
however, that the uncertainty in the signal could be significantly 
larger than the formal error because of the complexity in deter- 
mining the background (worst case 50 %). 

In |Grimes et al.|f2007| ), the same data set was reduced us- 
ing CalFUSE v3.1.8. The authors here chose to treat Haro 11 
as a point source, after they were certain that all the signal was 
included in the smaller aperture. The method to use the point 
source option can be well motivated also for galaxies, since this 
works to improve the S/N in two ways. Not only are the data 
corrected for astigmatism, but the smaller extraction region also 
minimizes the error introduced by the background model. To 
reduce the contribution from airglow lines and solar scattered 
light, only the night portion of the data was used, reducing the 
exposure time from 16 ks to 12 ks. The IDF-files from the differ- 
ent exposures were cross-correlated with each other and coadded 
before the background fit and extraction of the ID spectrum. No 
LyC signal could be seen by visual inspection for either 2D de- 
tector segment. The extracted ID spectrum gave a negative flux 
in the LyC for the SiC IB spectrum, but a weak positive flux in 
the SiC 2A (see Grimes et al. 2007, Fig. 7). The authors derived 
an upper limit to the escape fraction of f esc <2%. 



3.1.3. Analysis of CalFUSE v3.2 

Since the two previous articles yielded such different results, 
and considering the known problem of CalFUSE often over- 
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estimating the background at short wavelengths for faint targets, 
a detailed analysis of the Haro 1 1 data set reduced with CalFUSE 
v3.2 was undertaken. The analysis is not a comparison of the two 
previous articles, but instead describes the performance of the 
final version of CalFUSE when applied to low S/N data using 
different settings such as point source, extended source, normal 
extraction , an d optimal extraction. The analysis is presented in 



Appendix ^ 



We also discuss whether Haro 1 1 should be treated 
as an extended or point source, and take a look at the special case 
where the signal might originate from off-center regions. 

Since the analysis is quite technical, we only present the 
main conclusions here and refer to Appendix |A| for details. 
i) We cannot rule out that Haro 1 1 might exceed the limits de- 
fined for an object to be treated as a point source. We know that 
star-forming regions exists out to large distances, even outside 
the LWRS aperture (Fig. [2}. In addition, the spectral height of 
the Haro 1 1 continuum is found to be quite extended, roughly 
half that of the filled aperture airglow emission. From Fig. A.l 



it seems as if some of the signal actually might fall outside the 
point source aperture. 

if) Using the point source option, the data are corrected for astig- 
matism in the optical design. This is made using a model based 
on a true point source and is applied to the brightest source 
within the aperture. When applied to a galaxy where the emis- 
sion might come from several regions of varying intensity with 
wavelength, we have no way of knowing if this is done in a cor- 
rect way. For faint targets, the correction to the SiClB spectrum 
could possibly even be applied to the slope of the background 
itself, providing the brightest pixels in parts of the spectrum. 
Hi) Using the point source option, the extraction of the ID spec- 
trum is by default made with the optimal extraction routine, 
meaning that central pixels are assigned a higher weight. Since 
we cannot rule out that any LyC emission might originate from 
off-center regions, such a treatment needs to be avoided and all 
pixels should be given equal weight. 

iv) Non-simultaneous estimates of the background are insuffi- 
cient for faint targets. The template background files are seen to 
vary in an uncorrelated way with the Haro 11 data (Fig A.2i, 
which can cause large uncertainties in the final measured flux. 
The effect is most severe at the shortest wavelengths, often re- 
sulting in an overestimated background and negative fluxes. 

v) Even though the extended source option was found to be tech- 
nically the more reliable option, the larger extraction region ap- 
plied with this method also increases the uncertainty in the back- 
ground model. 

The above conclusions found from the analysis in 
Appendix[A] led us to develop a new background routine for re- 
ducing the FUSE data of faint targets, and to include both point 
source and extended source reduced data in our analysis. 

3.2. A new background model 

CalFUSE was not originally designed for dealing with objects 
with very low photon counts (Dixon 2007, private communica- 
tion). The main problem lies in how the background is modeled, 
particularly in the SiC 1A and IB detector segments where the 
spectra lie on the steeply rising background at the lower edge 
(see Fig. A.l and Fig. 4.1 in DH). In Appendix [A] we show that 
non-simultaneous background estimates with a predefined shape 
are not applicable when dealing with low signal-to-noise data. 
The approach here is instead to use the more standard method of 
directly modeling the background from the data, which also has 




The Appendix is only available in the online journal. 



Fig. 3. The IB (left) and 2A (right) detectors plotted in the cross- 
dispersion direction. The data were reduced as extended source, 
and were binned with a factor 128 in x and 8 in y, to a 128x128 
pixel sized image. The data are summed over the pixel inter- 
vals x=10-20 (upper), x=50-60 (middle) and x= 100-1 10 (lower). 
The original signal is plotted in gray, with three different back- 
ground fits overplotted in black: A X =A V =0.5 (solid), A x =A y =2.0 
(hatch/dot, and A A =A V =6.0 (hatched). The SiC spectra are cen- 
tered at y pixel ~10 for IB and y pixel ~ 50 for 2 A, while the 
LiF spectra are centered at ~ 70 for IB and ~ 92 for 2A. The 
narrow dips and peaks that are the most prominent on the IB 
detector are artifacts of pipeline processing (Dixon 2007, private 
communication). 



the advantage of dealing with features from stray and scattered 
light described in DH, Section 7.2.1. 

The new background model is constructed as follows. First 
we run CalFUSE v3.2 using both the extended source and point 
source options to obtain the IDF files, and combine the different 
exposures as recommended in Dixon et al. ( 2007 ). Since all sub- 
images are co-added before the fit, we automatically deal with 
any time variations in the data. All the photon events are then 
reconstructed using the final pixel x- and y coordinates from the 
IDF, and the data are binned to a 128x128 pixel sized image. 
The binning is made under the assumption of background conti- 
nuity and smoothness, and allow a higher signal when creating 
the background model. A mask of matching size is then pro- 
duced that excludes the regions where the signal lies and the 
edges of the detectors. Since it is difficult to achieve a good fit 
to the steep edges at the upper and lower parts of the detectors, 
and the drop-off outside these, we mask out everything outside 
the maximum pixel value in the edges. For segment IB, the re- 
maining pixels between the SiC aperture and the edge are too 
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Fig. 4. The signal in the 2A detector segment plotted over the 
whole x range, summed over row pixels y=37-43 (left panel) 
and y=97-103 (right panel). The left panel samples the back- 
ground and background fits just below the SiC spectrum, while 
the right panel samples the background and background fits just 
above the LiF spectrum. The original signal is plotted in gray, 
and three different background fits are overplotted: A x =Aj,=0.5 
(solid), A A =A V =2.0 (hatch/dot), and A. V =A V =6.0 (hatched). 



few in number (sometimes even zero) to be able to perform the 
fit. Therefore, the point source extraction window is imported 
and used for defining the mask also for the extended source 
spectrum (but not for the extraction of the ID spectrum). The 
background is then modeled by fitting a smooth surface to the 
128x128 masked image using optimal filtering techniques. We 
achieve this by solving a regularised optimisation problem 



J] M X ,y [D X ,y - F X ,yf + A, [F X+h y - F X ,yf 

x,y x,y 

+A.v ^ \Fx,y+\ - F x , y f = min 



where M x y , D xy , and F x>y are the mask, the data, and the model 
value at image "pixel" x, y. The constants A v and A v control the 
smoothness of the model in the x and y directions, and are se- 
lected by the user. The model is computed by solving a sparse 
system of linear equations with a band-diagonal matrix. The 
solver is implemented as an IDL routine opt Jilter J.d.pro writ- 
ten by N. Piskunov, which is part of the Reduce package^* The 
Reduce package is described in Piskunov & Valenti (2002), but 
since the opt jilter J.d.pro was developed later, a description of 
the routine was not included in this paper. 

In Fig. [3] the background fit is shown at three different x- 
positions for both detector IB (left panels) and 2A (right pan- 
els), each with three different choices of smoothing factors. The 
combination of A A , A v that yielded the best fit to the background 
with avoidance of the signal areas was found to be 0.5, 0.5 for 
detector segment IB, and 2.0, 2.0 for 2A. The lower smoothing 
factor for IB was necessary in order to try and fit the steep slope 
of the background in the SiC spectrum, but even using this loose 
constraint the fit failed to reproduce the shape to a satisfactory 
level. The reason for this is simply that after the mask is applied 
there are not enough data points left on the slope to perform the 
fit, even though the more narrow point source extraction win- 
dow was used. For the 2A segment on the other hand, the SiC 
spectrum falls on a region of the detector where the background 
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4 The Reduce package and opt_filterJ2d.pro routine can be found at: 
|http://www.astro.uu.se/~pis kunov/RESEARCH/REDUCE/r educe.htrnl| 
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Fig. 5. The pulse height distributions for three different galax- 
ies observed at different voltage settings throughout the mission. 
The pulse height distribution of the 920-980 A wavelength re- 
gion in the SiC 2A spectrum is plotted in gray, with the back- 
ground distributions from an equally large area under the spec- 
trum plotted in black. Upper: Tol 0440-381, Middle: Haro 11, 
Lower: M 100. The drift between the two distributions show 
that the differential loss of sensitivity for the spectrum and back- 
ground regions could have a significant effect for M 100, which 
was observed late in the mission after several detector voltage 
changes. In addition, for Haro 11 and Tol 0440-381 there is a 
slight shift, indicating that there might be some effect present. 



is much more well behaved. The fit was found to be rather in- 
sensitive to the choice of smoothing parameter in the SiC spec- 
trum, while in the LiF spectrum closer to the upper edge it is 
more dependent. In Fig. [4j this is illustrated for the whole x- 
range for y-pixels=37:43 (just below the 2A SiC spectrum) and 
y-pixel=97:103 (just above the LiF 2A spectrum). The flux was 
calibrated using an interpolation between the two CalFUSE ef- 
fective area files closest in time to the Haro 1 1 observations. The 
error in the background fit for the SiC 2A spectrum is shown as 
the hatched lines in Fig. [8] and the error is propagated through 
to the final ID spectrum and the S/N calculations. 

Since the SiC IB fit was found not to be sufficiently good, a 
different approach was tried for this segment. First we reduced 
the binning of the IDF data, and produced a 256x256 pixel im- 
age. This made it easier to locate the separation between the 
slope of the background edge and the signal in the SiC LWRS 
spectrum. Instead of masking out the point source region as be- 
fore, we instead fit the background slope under the signal with a 
Gaussian. In a region of four pixels over the whole x range, eight 
points on the slope were interactively recorded to which a fourth- 
order polynomial fit was applied. After visual confirmation of its 
closeness, a given fit solution was saved to the image from which 
the background was modeled using the opt_filterJ2d.pro routine 
as before (A A =A V =2.0). To minimize the impact of the user, nine 
of these backgrounds were produced from which a mean image 
was made. The SiC IB extended source 2D image summed over 
the LyC region with the fit overplotted is shown in the lower 
panel in Fig. [6] The ID spectrum with the mean background 
overplotted can be seen in Fig. [9] Here, also the standard error of 
the background fit calculated from the nine images is displayed 
in lighter gray. No similar attempt to fit the background of SiC 
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IB with point source reduced data was tried since the separation 
of the background slope and the signal is very difficult to iden- 
tify, and the shape of the slope of the background itself is also 
hard to estimate (see Fig. |A. 1) >. 

An effect that is not taken into account in the construction 
of this background model is the gain sag of the detectors, which 
might lead to an underestimate of the flux. Owing to the loss 
of sensitivity with time and exposure, the detector voltage was 
increased several times during the mission resulting in a shift 
towards higher pulse height values (DH, Table 7.4). The il- 
luminated parts of the microchannel plates lose more of their 
efficiency over time than the unilluminated parts (IH, Section 
4.4.2.1). Since we mask out the extraction regions and fit the 
background to the unilluminated parts, this means that the back- 
ground in our model may be somewhat overestimated. The Haro 
1 1 data were obtained early in the mission, when this effect was 
not expected to be significant. To investigate the effect, we plot- 
ted the pulse height distributions of three galaxies observed with 
three different detector voltage settings. Tol 0440-381 (2000-12- 
13) was observed with the initial voltage setting, Haro 1 1 (2001- 
10-12) after the second voltage change, and M 100 (2004-02-18) 
late in the mission when the voltage had been changed several 
times. In Fig. [5] the pulse height distributions for both the 920- 
980 A wavelength region in the SiC 2A spectrum (gray) and the 
background in an equally large area below the spectrum (black) 
is shown. For M 100, the shift is significant between the two dis- 
tributions, indicating that the differential loss of sensitivity can 
have a significant effect for this target. For Haro 1 1 and Tol 0440- 
381, there is a slight shift between the two distributions, and the 
differential sensitivity could also have some effect for these tar- 
gets. If an effect is present, it will only work in the direction of 
more conservative LyC flux measurements. 




100 150 
y pixel 



4. Results 

4.1. The Lyman continuum measurements 

In Fig. [6] the night-only signal has been summed over the LyC 
range (920-929 A) in the SiC spectrum for the 2A (upper panel) 
and IB (lower panel) IDFs, respectively. A weak excess can be 
distinguished in both spectra above the background fit plotted 
in black. In an effort to evaluate this signal and constrain the 
contributions from possible noise sources, we took a closer look 
at the data concerning airglow and solar scattered light, count 
rates, and Milky Way absorption lines. 

4.1 .1 . Lines affecting the Haro 1 1 LyC region 

Terrestrial airglow and solar scattered light emission unfortu- 
nately affects the wavelength region where we wish to mea- 
sure the LyC of Haro 11. The airglow emission was studied 
by |Feldman"et~aT1 ( p00T] i, using FUSE. In Fig. 7 of that article, 
they present spectra obtained at three different limb angles (limb, 
disk, and dayglow spectrum). The limb and disk exhibit distinct 
spectral features that are of the order of 30 times stronger than 
that of the dayglow spectrum. If these features are present in a 
spectrum it would be obvious, but neither can be identified in 
the Haro 1 1 data. Some dayglow features (airglow) are however 
clearly present at longer wavelengths (e.g. OI /1989A and Ly/3 
A 1025 A), and in some cases in both the orbital day and night 
data. Although the airglow lines are weak or even invisible in 
the LyC wavelength region, we must consider a possible contri- 
bution from these. The airglow emission lies on top of the Milky 
Way Lyman absorption lines, so the emission and absorption are 



Fig. 6. The 2 A (upper) and IB (lower) segments plotted with 
night-only data. The signal summed over the LyC range (920- 
929 A) in the SiC spectrum is plotted in gray, with the new 
background fit overplotted in black. The vertical hatched lines 
indicate the extraction regions for the LWRS extended source 
apertures. The SiC 2 A spectrum is centered at y pixel ~50, while 
the SiC IB is centered at pixel ~20. A weak excess can be seen 
in the LyC region for both spectra. 



seen superimposed. The affected Lyman line regions have a lim- 
ited extent in wavelength, so their contribution can easily be ex- 
cluded when measuring the LyC of Haro 1 1 . 

The solar scattered light is caused by solar light scattered 
within the SiC channel baffles (DH, Ch 7.1). This emission can 
be avoided by using night data only, so the data used for the 
analysis in this paper should be unaffected by this scattering. 
By inspecting the day and night data separately, we did however 
locate a residual emission just below 930 A, which is likely a O I 
airglow line. This region was excluded together with the other 
airglow regions. All the excluded regions are displayed as gray 
vertical bars in Fig. [8] 

4.1.2. Count rate versus limb angle 

The count rate columns in the IDF files (extension 3) are pro- 
duced by excluding the airglow and scattered light regions de- 
fined by the CalFUSE airglow calibration files. To check the per- 
formance of this filtering, a test was performed on the Haro 11 
data. Following |Grimes et al.| (see|2007 Fig. 12), we plotted the 
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Fig. 7. The count rate for SiC 2A plotted as a function of Earth 
limb angle for three different wavelength windows. Airglow 
photons have been excluded. Orbital day data are plotted in the 
left panel and orbital night data in the right panel. The plots do 
not show any obvious trend for increasing count rate at lower 
limb angles neither for the day or night data, indicating that the 
pipeline is successful in sorting out airglow and scattered light 
photons. The horizontal lines show the background count rate for 
the 920-930 A window. A faint excess in the LyC signal above 
the background count rate can be seen for most limb angles. 



count rate as a function of limb angle for SiC 2A. The idea is that 
if the airglow and scattered light photons are successfully sorted 
out, the count rate will be constant with limb angle. Another in- 
teresting possibility for this plot is to see whether there is any 
excess in the Lyman continuum. The count rate for the 920-930 
A wavelength region would in that case have to lie above the 
background count rate. 



Grimes et al. ( 2007 ) plotted the count rate for Haro 1 1 in 



three wavelength windows, 920-930 A (LyC region), 930-940 A, 
and 940-950 A, and found a clear trend for a decreasing count 
rate towards higher limb angles in all three bins. They also ar- 
gue that the 920-930 A bin for high angles is close to the nom- 
inal intrinsic background count rate, implying no LyC leakage. 
It was unclear from their discussion why in their plot the other 
two wavelength windows are also close to the intrinsic value at 
higher angles, which is remarkable at least for the 940-950 A bin 
where there is a strong stellar continuum flux from Haro 1 1 (see 
Fig.|8j. 

In Fig. [7] we reproduced this plot for SiC 2A. Since all the 
Haro 1 1 orbital day data were obtained for limb angles <38 de- 
grees, and all the night data were obtained for limb angles larger 
than 38 degrees, we plot the count rates for these components 
in two different panels. As expected, there is a clear distinction 
between the data for the two components, there being a higher 
count rate overall during orbital day. From the plot, we can draw 
the conclusion that, if we exclude the noisy data point at the low- 
est limb angle, there is no clear trend for higher counts at lower 
limb angles for either the day or night data. This suggests that the 
airglow regions indeed are successfully filtered out by the cali- 
bration files of the pipeline, if that option is set by the user. The 
difference between day and night data can likely be attributed to 
a general increase in the background due to scattered light. 

We also wish to explore the possibility of drawing conclu- 
sions about possible LyC leakage from this plot. The Haro 1 1 ob- 



servations were obtained under different conditions from when 
the nominal intrinsic background was measured. The latter was 
derived at the end of the mission in 2007 with the telescope 
doors closed, and was averaged over the whole detector (IH, 
Section 4.4.3.1). This means that the two background measure- 
ments should differ in several ways, such as stray light (which 
was not present when the doors were closed), the higher voltage 
on the detectors at the end of the mission, and the solar cycle be- 
ing different (Haro 1 1 observations were taken well before solar 
minimum, while the nominal background was measured during 
a solar minimum). Instead of plotting the nominal background, 
we therefore extracted the background derived in the modeling 
of the Haro 11 data (Section 3.2 1. The background count rate 
for the 920-930 A window is plotted as a horizontal hatched line 
in Fig. [7] It is not a function of limb angle, but an averaged value 
in the 920-930 A window scaled by the exposure time and size 
of the region excluded due to airglow. A weak excess in the LyC 
above the background can be seen for most limb angles in the 
night panel. The day data is only 25 % of the total observation, 
and the larger errors make it difficult to draw any conclusion 
from this panel. 

We also note that the 940-950 A bin as expected has values 
high above the background, which in this wavelength window is 
0.070 (day) and 0.047 (night) cnts s" 1 . 



4.1 .3. The 1 D spectrum 
The SiC 2A spectrum 

The ID SiC 2A spectrum of Haro 1 1 using the new background 
model can be seen in Fig. [8] In the two upper panels, the data 
were reduced as an extended source (EC), while in the two lower 
panels the data were reduced as a point source (PC). The left 
spectra show the total signal (black), the background (dark gray), 
and the error in the background fit (light gray). The right spectra 
show a close-up of the Lyman continuum region, with the back- 
ground subtracted from the total signal (i.e. the final flux density 
of Haro 11). Since day data are affected by airglow and solar 
scattered light emission from Oi and Hi (DH, Ch 7.1), only the 
night portion of the data was used. The night side Hi airglow 
regions still need to be masked out because of the presence of 
interstellar absorption possibly blended with geocoronal airglow 
emission, and are marked as gray vertical bars. 

The LyC flux density (/900) was measured on the airglow 
free regions below the Haro 1 1 Ly -limit at 930.3 A, which gave 
a total wavelength interval of 7 A. For the EC option, the flux 
density was found to be /9oo,ec=4.2 x 10~ 15 erg s _1 cirT 2 A -1 
with S/N=3.8. For the PC option, the flux density was found to 



be /9oo,/>c=2.8 x 10 erg s 



with S/N=2.7. 



The Milky Way extinction at 920 A (corresponding to Haro 
11 rest frame 900 A) along the line of sight is A92o,mw=0.1 17 
(Bergvall et al. 20061, and the corrected flux densities are 
/900,£C=4.6 x 10 ' 15 and / 9 oo,pc=3.1 x lO" 15 erg s" 1 cirr 2 A" 1 . 

The errors were calculated as follows. To evaluate the error 
in the background fit, regions above and below the SiC 2A 
spectrum on the 2D IDF were defined with total area equaling 
that of the aperture. For each resolution element the standard 
errors are calculated for these two regions, with respect to the 
original data and the background fit. The total standard error 
for each airglow free window is then calculated from both the 
error in the background fit and the photon statistics (Gaussian) 
in the ID spectrum. Finally, the weighted mean of the airglow 
free windows is calculated for both the signal and the error in 
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Fig. 8. The Haro 1 1 SiC 2A spectrum. The upper panels show the data reduced as extended source (EC) and the lower panels as 
point source (PC). The binning is 0.23 A. Left: The total signal is plotted in black, the new background in dark gray and the error in 
the background fit in light gray. We also indicate some intrinsic Haro 1 1 Lyman lines, of which only Lye at 957 A falls in a region 
uneffected by geocoronal or solar scattered light emission. Right: A closer view of the LyC region. Here, the background has been 
subtracted from the total signal, and the plot hence show the final flux density of Haro 1 1 . The airglow regions are marked with gray 
vertical bars. 
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Fig. 9. The Haro 11 SiC IB spectrum for the extended source option. The binning is 0.26 A. Left: The total signal is plotted in 
black, the new background in dark gray and the standard error in the background in light gray, we also indicate the intrinsic Haro 
1 1 Lyman lines down to the Lyman limit. Right: The final flux in SiC IB overplotted with the final flux in SiC 2A (gray). We label 
the Milky Way Lyman absorption lines, for which some geocoronal emission also can be seen at the center. 



the LyC region, where the weight is calculated from the size of 
the window and the effective area. In the S/N calculations the 
unidentified but correlated absorption features between EC and 
PC seen in the right panels of Fig.|8]were treated as noise. 
The SiC IB spectrum 

In Fig. [9] the ID spectrum of SiC IB for the extended source 
option is shown. In the left panel, the total flux is plotted in 
black, the background in gray, and the standard error in the 
background in lighter gray. There is a significant difference 
between the background for SiC IB and SiC 2A. While the 
2A background varies little and smoothly over the whole 
wavelength interval, the IB the background varies dramatically. 
In the right panel, the final flux (total flux minus background) 



is plotted in black for SiC IB and, for comparison, the final 
flux in SiC 2A is overplotted in gray. The spectra seems to 
agree fairly well, but the flux is slightly stronger over all in 
2A. It is especially interesting to see is that below ~ 920 A 
the flux stabilizes around zero, which might be expected since 
the Milky Way Lyman forest is thick at these wavelengths. An 
excess in the LyC is also noted, with /9oo,ib=3.2 x 10~ 15 erg s~ ! 
cm -2 A -1 with S/N=2.6. Owing to the complex nature of the 
background surrounding the SiC IB spectrum, the method for 
calculating the errors in the background fit differs from that used 
to calculate the errors in the SiC 2A fit. Here the errors in the 
fit are calculated as the standard error in the nine backgrounds 
described in Section 13.21 
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Including day data 

As a comparison, we also looked at the day+night data in the 
SiC 2A. The corrected fluxes were found to be /90o,£C=5.1 x 



10 erg s 



with S/N=3.5, and / 90 o„pc=3.1 x 10" 



erg s cm A with S/N=2.7. Both values are within one cr 
of the night-only results, indicating that scattered light during 
the orbital day has a minor influence on the data. 

The final flux density 

In the calculations of the escape fraction presented in the next 
section, we adopted the combined flux density from SiC IB 
and 2A extended source for night data only. The compatibility 
between the two spectra was evaluated on the 920-965 A 
region, avoiding the airglow features around 930 A. The null 
hypothesis that the two spectra are statistically comparable is 
supported at a confidence level of more than 90 % according 
to the Kolmogorov-Smirnov test, and a reduced chi-square test 
of the difference between the two spectra gave a result close to 
unity, X 2 red = l - 07. The final LyC flux is / 9 oo=4.0 x 10~ 15 erg s" 1 
cm" 2 A" 1 (S/N=4.6). 



4.1 .4. Zero flux level 

To validiate our results, we need to check the zero flux level. The 
narrow Milky Way absorption line C m A911A should be satu- 
rated at line center, if the line is resolved and free of solar scat- 
tered light emission (which should be absent when night-only 
data are used). In Fig[l0j where the SiC 2A night-only spectrum 
has been zoomed in on the line, it does seem as if the flux zero 
level is accurately calibrated using our background model. The 
intrinsic Haro 11 hyd redshifted to 969.2 A, also appears to be 
saturated at line center. 

Another possible way to calibrate the zero level was pre- 
sented in |Zheng et al.] pOlM] ), where the faint QSO HE 2347- 
4342 was analyzed. Using night-data only, saturated intergalac- 
tic He ii lines were fitted with a smoothly varying polynomial 
to determine the zero level as a function of wavelength. For 
the Haro 1 1 data set however, this approach was not an option 
since the Ly-lines are also affected by airglow when night-only 
data are analyzed. In Fig 10 both the airglow emission and the 
broader interstellar absorption can be seen for Lyy at /1972.5A. 
The airglow emission is observed down to ~930 A, and below 
that no clear lines can be observed indicating that the airglow 
emission and interstellar absorption are well blended and some- 
times even cancel each other. 



4.2. The LyC escape fraction 

Following [Deharveng et al. (2001), the absolute escape fraction 
f esc is defined as the ratio of the observed flux density of 900 
A photons escaping the galaxy without being absorbed, fgoo, bs> 
to the 900 A flux density actually produced by the stars in the 
galaxy, / 900 Jnt •> 

given by 



fesc — 



A 



900,obs 



Aooj 



(3) 



To estimate the absolute escape fraction of LyC photons, we 
need to know the total number of ionizing photons produced 
in the galaxy. A detailed way of deriving the intr insically pro- 
duced 900 A flux for Haro 1 1 was presented in |Hayes et al. 
( |2007| l. With high resolution HST images, they were able to 
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Fig. 10. Part of the Haro 1 1 SiC 2A extended source spectrum 
plotted with binning 0.08 A and night data only. The Milky Way 
C m A977 A line is seen to be saturated, indicating that the zero 
flux level is accurately calibrated using the background model 
presented in this paper. In addition, the Haro 1 1 Ly<5 line red- 
shifted to 969.2 A seems to be saturated at the line center. The 
Milky Way absorption Ly-lines are however affected with air- 
glow emission, and cannot be used as a calibration tool for Haro 
1 1 (here Lyy at ^972.5A). 

model the fgoojm at eac h binned pixel based on the normaliza- 
tion of the best- fit population synthesis models. Integrating the 
ionizing flux over the whole galaxy, they derived a total flux den- 
sity of /9oo.hh-12.3 x 10~ 14 erg s cm~ 2 A -1 . Using this value, 
we derive an escape fraction for Haro 1 1 of / csc =3. 3+0.7 %. 



4.3. The C 11 A1036A escape fraction 



|Heckman et al. ( 2001 1 presented a method for indirectly measur- 
ing the LyC escape from galaxies. The method is based on the 
expectation of residual emission in the center of the C 11 A 1036 A 
interstellar absorption line if the target is a LyC leaker. Five lo- 
cal FUSE galaxies were observed in this pilot study, and they 
were all found to be essentially black in the Cn /11036A line. 
Grim es et al.| ( |2009| l and |Heckman et al.| ( |201 1[ ) increased the 
sample size with a factor ^3 using FUSE archival data. Neither 
found strong evidence of any large leakage, although three of the 
galaxies were found to have rather high upper limits, / esc <10- 
18%. 

Haro 1 1 was also included in the sample. The ratio of the line 
core flux of C 11 /U036A to that of the continuum was found to 
be Rcv/=8 %. Applying CalFUSE to detector segment 1A with 
the background model presented in this paper, we measured a 
similar ratio of Rc//=8.5 %. After correction for the far-UV ex- 
tinction, they arrived at an upper limit of f esc ,LyC < 0.6%. At first, 
this may seem to contradict our results of f esc , but we claim that 
this is not necessarily the case. 

Since it has been shown that Haro 1 1 is optically thick at 900 
A if the absorbing gas is diffusely distributed, any emitted LyC 
photons must escape following the picket fence model. Because 
of this, there are two effects that must be considered. Firstly, 
there are two components in the spectrum, one that is attenu- 
ated by dust and one that is not. Differential extinction between 
the unattenuated Lyman continuum at 900 A and the attenuated 



spectrum at the wavelength used by |Grimes et al. (2009 1 to cal- 
ibrate the UV flux, i.e. A=l 150 A, plays an important role. The 
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escape fraction can only be the same at 1036 A as at 900 A if the 
spectral shape of the reddened part of the spectrum is the same 
as that for the part leaking through the fence, which is an un- 
likely scenario. Following the SMC extinction curve published 
by |Cartledge eTaT^ ( [2005) 1 and A 960 =2.9 m ( [Bergvall et al.|2006| , 
we obtain a difference in extinction between 1036 A and 900 A 
corresponding to a factor of 1.3. Secondly, we are probably de- 
tecting two different stellar populations, one whose light is ob- 
served through the holes and is dominated by very young stars 
(1-3 Myr), the other whose light originates from the surrounding 
regions which contains a more mixed population. The far-UV 
slopes of the two populations are therefore likely to vary, and 
the fluxes deviate significantly when extrapolating to LyC wave- 
lengths. Since most of the light of Haro 1 1 detected in the far-UV 
is enclosed within the LWRS aperture, the mean age of the stars 
in the UV spectrum will be significantly higher than 1-3 Myr. 
Berg vall et al.| ( [2006| derived a best-fit age of 14 Myr using the 
models of synthetic spectra for star-forming galaxies by [Robert] 
|et al.| ( |2003[ ). This must be considered as a lower limit taken into 
account that the leaking population has not been subtracted prior 
to the fit, and that the models only include spectra for stars with 
ages <15 Myr. In their study of young stellar clusters in Haro 1 1, 



Adamo et al. ( 20 1 0)> claim that the burst has lasted for 40 Myr. 



Two plausible models demonstrating the difference in spec- 
tral energy distribution between /11036A and the Ly-limit with 
both these effects taken into account, have been summarized in 
Fig. 11 In the upper panel, the Starburst99 model (Leitherer 



et al.|199"9[ ) is plotted with the leaking population (3 Myr, con- 
stant star-formation rate (SFR)) in solid line, and two non- 
leaking populations (hatched: 20 Myr with constant SFR, dotted: 
20 Myr, instantaneous burst). The lower panel obtained with the 



Zackrisson et al. ( 200 1 )> model, essentially show the same thing 



but the non-leaking model is plotted for a 20 Myr population 
with a constantly declining SFR over a timescale of 1 Gyr. The 
models are for stellar populations with 20 % solar metallicity, 
similar to that of Haro 1 1 . All non-leaking models have been 
corrected with the same extinction as observed in Haro 1 1 . The 
ratio of the leaking to non-leaking populations is most sensitive 
to the extinction and the age of the leaking population. While 
younger stars would lead to a higher correction factor, a domi- 
nating older population would make the effect smaller. Raising 
the age of the non-leaking population to 50 Myr has however 
almost no effect at all. 
Both models in Fig. 



1 1 indicate that with the combined ef- 



fects, the extinction-corrected LyC flux density derived from the 
C n /U036A line should be corrected with a factor »4 from 1036 
A to 900 A. This value was derived using a reasonable model 
with a leaking population of age 3 Myr. Correcting the abso- 
lute escape fraction derived by Grimes et al. (2009), we arrive 
at fesc,LyC < 2.4%. Calculating with our measured ratio Rcv/=8.5 
%, and A%o=2.9 m , we arrive at the same value, 2.4 %. 



5. Summary and conclusions 

As of today, direct detection of Lyman continuum leakage has 
only been claimed for one galaxy in the local universe, Haro 1 1 
(Berg vall et al.|2006] l. In a similar analysis, using a later version 
of the FUSE data reduction pipeline, Grimes et al. (2007) could 
not confirm this detection. The aim of this work was to settle the 
question and perform a detailed analysis of Haro 1 1 data for the 
two channels covering the LyC region, SiC IB and SiC 2 A. 

In previous articles where z~3 galaxies were imaged in the 
LyC, it was shown that emitted ionizing photons can be found at 
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Fig. 11. Predicted spectral energy distributions in Haro 11 il- 
lustrating the picket-fence model. The stellar spectra are ob- 
tained from two models with 20 % solar metallicity - Starburst99 
( Leitherer et al.|1999 > (upper) and our model (Zackrisson et al.| 
|2001[ ) (lower). We assume that the leaking stellar population has 
a young age, 3 Myr, has had a constant SFR and is unaffected by 
extinction. The non-leaking component is assumed to be 20Myr. 
In the Starburst99 model, we look at two star-formation histo- 
ries of this component, an instantaneous burst and a model with 
continuous star-formation. In the Zackrisson et al. model, we as- 
sume an exponentially declining SFR on a timescale of 1 Gyr. 
The old component has been reddened by assuming the same 
extinction as observed in Haro 11. The Zackrisson et al. model 
spectra have been smoothed with a triangle smoothing to reduce 
the differences in spectral resolution between the stellar spec- 
trum (low resolution) and the nebular spectrum (high resolution). 



some distance from the center of the galaxies ( |Iwata et al.|2 009 ; 
[Nestor et al.||201l) . While there is no possible way of acquir- 
ing LyC images of Haro 11, we do know that this local galaxy 
has star-forming regions out to, and beyond, the 30x30"LWRS 
aperture of FUSE. Even though the spectral height of Haro 1 1 at 
the wavelengths of the non-ionizing continuum was found to be 
close to the size of a point source target, a simple model indicates 
that the LyC emission might be seen as extended by FUSE. The 
analysis in this paper was therefore made with both the extended 
source and point source pipeline options. 

The detailed analysis of the reduction pipeline in 
Appendix [A] shows that one has to be very careful using the 
point source option of CalFUSE v3.2 for galaxies and poten- 
tial Lyman continuum leakers where the signal might originate 
from off-center regions. We discuss the importance of carefully 
choosing the cross-dispersion centroid for the extraction regions, 
the risk of astigmatism corrections, and the skewed weighting of 
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central pixels as the ID spectra are extracted using this option. 
The extended source option on the other hand, although found to 
technically be the more reliable option, increases the uncertainty 
introduced by the background fit. 

The non-simultaneous background model used by the reduc- 
tion pipeline was found to be inadequate for the low S/N spec- 
trum of Haro 1 1 . The shape of the template background file can 
be seen to vary in a way sometimes uncorrected to the Haro 1 1 
data, and the effect is most severe at the shortest wavelengths. A 
new method for subtracting the background was therefore devel- 
oped, where the background is modeled by fitting a smooth sur- 
face directly to the detector using optimal filtering techniques. 
The new background model was found to be very accurate for 
the SiC 2 A spectrum, but not for the SiC IB spectrum positioned 
on the steep slope of the detector background. For this spectrum, 
the model had to be modified by fitting a Gaussian to the slope of 
the background of the target signal. No attempt to fit the SiC IB 
point source reduced data was made, owing to the difficulty in 
distinguishing between the background slope and target signal. 
A weak LyC excess in the 2D data above the background could 
be identified in the night-only data of both SiC IB and SiC 2A. 

Plotting the count rate versus limb angle displayed no trend, 
indicating that the calibration files of the pipeline successfully 
filters out airglow and solar scattered light photons if this option 
is set by the user. In addition, by comparing the count rate in the 
LyC range (920-930 A) to the corresponding background count 
rate, we found an excess in the LyC for most limb angles. 

Applying CalFUSE v3.2 to the Haro 11 data set in addition 
to the new background model, we extracted the ID spectra. The 
LyC signal observed in the 2D data and count rate plot could 
also be confirmed here. Using the extended source SiC IB and 
SiC 2A night only data, the flux density in the LyC between the 
airglow regions was found to be /9oo=4.0 x 10~ 15 (S/N=4.6) erg 
s _1 cirT 2 A -1 . The errors were calculated including both photon 
statistics and the uncertainties in the background fit. From this, 
we derive an absolute escape fraction of / esc =3.3±0.7 %. 

Since it has been shown that Haro 11 is optically thick at 
900 A if the absorbing gas is diffusely distributed, any emit- 
ted LyC photons must escape following the picket fence model. 
Previous works where the LyC flux was estimated from a resid- 
ual emission in the centre of the C n /U036A interstellar absorp- 
tion line, failed to take into account two important effects using 
this model. These effects of differential extinction and the differ- 
ent spectral energy distributions of the leaking and non-leaking 
populations, should in the case of Haro 1 1 cause a factor of four 
difference in flux between 1036A and the Lyman limit. Thus we 
derive an escape fraction of < 2.4 % using the C n /11036A line, 
in agreement within the errors with our measurement obtained 
directly from the Lyman continuum. The correction factor (and 
fesc) will be larger if the age of the leaking population is younger 
than the assumed 3 Myr. 

Thus Haro 1 1 continues to be the only galaxy in the local uni- 
verse for which escaping LyC photons have been detected. The 
low escape fraction is not unexpected by model predictions. One 
important parameter in these models is the mass, where most 
models seem to agree on a smaller f esc for galaxies with higher 
masses such as that of Haro 1 1. In addition, the increase in metal - 
licity with redshift and the lack of PopIII stars greatly reduce the 
probability of finding large escape fractions from local galaxies. 
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Appendix A: Using CalFUSE v3.2 

Performing LyC measurements requires data acquired while 
working at the very limit of instrument sensitivity. We have 
therefore dedicated this appendix to a detailed analysis of the 
performance of the reduction pipeline on the Haro 1 1 data set. 
Haro 11 is one of several faint FUSE targets for which the 
CalFUSE background subtraction has resulted in negative flux 
at short wavelengths, and the results found here may be of inter- 
est to other similar investigations. We also approach the special 
case where the emission might originate from extended or off- 
center regions of the targets. 



A. 1. Haro 11, extended or point source? 

This is a fundamental question when the reduction approach is 
to be decided. From the ACS/SBC UV image in Fig. [T] we see 
that the separation between the two strongest UV emitting re- 
gions is ~5". In the header of the Haro 1 1 FUSE data, the key- 
word APER_PA=99.4, which means that the distance between 
the two regions are almost maximized in the cross-dispersion 
direction. The LWRS aperture of FUSE is 30"x30". Only if 
the image height is roughly half the aperture size or more, will 
there be a significant effect on the height of the spectral image. 
Therefore, when only considering the far-UV ACS/SBC image, 
should Haro 1 1 pass as a point source. However, Fig.[TJcovers a 
square 15", and we know that star-forming regions do exist out- 
side this from Fig. [2] In the Chandra image of Haye s et aL] (see 
|2007j , Fig. 6) with the far-UV contours overp lotted, the X-ray 
emission is seen at least out to a diameter of 25". 

To check whether Haro 1 1 is consistent with a point source 
in the spectral image of the LWRS spectrum, we performed the 
following test. The astigmatic image height is smallest at 916 A 
in the SiC IB spectrum and at 926 A in SiC 2 A, that is, where we 
wish to measure the LyC. It is also close to reaching a minimum 
at Ly/3 in the LiF 1 A spectrum, where the continuum flux of Haro 
1 1 is strong enough to perform the test. The airglow Ly/3 fills the 
aperture, so if the Haro 1 1 continuum height is less than half of 
that it should pass as a point source. We find that the airglow Ly/3 
height is roughly 90 unbinned pixels, while the Haro 1 1 height 
just below and above Ly/3 is half of that, 45 pixels. We conclude 
that Haro 1 1 is a borderline case for being considered as a point 
source, and in our analysis we therefore considered both cases. 



A.2. Producing the IDFs 

Here, we reduced the Haro 1 1 data set with the last version of 
the pipeline, CalFUSE v3.2. The pipeline was run with both 
the point- and extended source options in order to compare the 
two methods. The processing trailer files were then scanned for 
warnings of which none were found. The IDFs and bad-pixel 
masks from the different exposures were coadded before the 
background fit and the extraction of the ID spectra was made, 
all according to the recommended strategy for low S/N spectra 
in |Dixon et al.| ( [2Q07] >. The IDFs were then cross-correlated be- 
fore they were combined (max correction 0.026 A). 
The cross-dispersion centroid 

The cross-dispersion centroids of the extraction windows were 
calculated using the cf_edit routine with option target events 
only, except for the SiC 1A and SiC IB where the contribution 
from the bright detector edge tends to move the centroids. For 
these two, the centroids were therefore determined with the op- 
tion user select after careful inspection. This treatment turned 
out to be especially important for the point source reduction. 
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Fig. A.l. Part of the SiC IB and SiC 2A IDFs seen in the cross- 
dispersion direction, summed over 920 - 950 A. The extraction 
windows for the LWRS aperture is plotted in vertical hatched 
lines. The cross-dispersion centroids were determined with user 
select for the IB files and target events only for the 2 A files, 
and the width shows the maximum width in the wavelength in- 
terval. In a) SiC IB reduced using the point source option, b) 
SiC IB reduced using the extended source option, in c) SiC 2A 
reduced as point source, and d) SiC 2A reduced as extended 
source. The point source spectra have been corrected for astig- 
matism. Comparing a) and b), it is clear that the effect of that 
is significant in the SiC IB spectrum. It is hard to determine 
whether the astigmatism correction has been made on the target 
data or on the slope of the background. The effect on the SiC 2A 
spectrum is not significant. The authors are not convinced that 
the flux from the target lies completely within the point source 
extraction windows for either IB or 2A. The dips and peaks to 
the right of the extraction regions in IB are artifacts of pipeline 
processing (Dixon 2007, private communication). 



In Fig. A.l the SiC IB and SiC 2 A IDFs were summed in the 
cross-dispersion direction over the wavelength region 920 - 950 
A. The extraction windows for the LWRS aperture are plotted 
as vertical hatched lines, with the width displaying the maxi- 
mum width within the wavelength interval. In a) and c) where 
the point source option has been used, part of the signal appears 
to fall outside the aperture. 
Weights and off-center emission 

In |Iwata et aL] ( |2009 1, LyC imaging was performed for a sam- 
ple of z^3.1 galaxies. It was found that the ionizing radiation 
in several cases was emitted from regions off-center from the 
non-ionizing emission. Looking at Haro 11, this possibility can- 
not be neglected since star-forming regions are also known to 
exist outside the HST/ACS far-UV image (see Sect. |2j, and the 
scenario would likely be applicable to other FUSE galaxies as 
well. The point source option that utilizes the optimal extrac- 
tion as standard routine must therefore be considered carefully 
when measuring possible LyC emission, not only because of the 
more narrow extraction window. This routine also extracts the 
ID spectrum using weights, where the central pixels in the ex- 
traction window receive larger weight than the ones closer to the 
edges. The extended source option on the other hand, is not as 
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Fig. A.2. The plots show the extracted Haro 11 total signal (galaxy + background -WEIGHTS) for SiC IB (upper row) and SiC 2A 
(lower row), over-plotted by the fitted background (BKGD) in dotted line. From left: a) and e) extended source setting, b) and f) 
point source setting with normal extraction, c) and g) point source setting with optimal extraction and WEIGHTS produced with 
background-subtraction routine turned off (BKGD imported from d) and h)), d) and h) point source setting with optimal extraction 
(corresponding to the treatment in |Grimes et al. ( 2007| )). In theory c) should be identical to d), and g) identical to h), but they are 
not as discussed in the text. The spectra are rebinned by 16 pixels=0.208 A. Only night data are plotted. 



sensitive to the way in which the centroid is chosen, and the ID 
extraction routine is not applied using weights. 
Astigmatism corrections 

When using the point source option, the data are corrected for 
astigmatism in the optical design. The model for this correc- 
tion is based on a true point source, and its applicability must 
be questioned when considering the possibility that LyC leakage 
might originate from off-center regions. The astigmatism cor- 
rection is performed on the location of the brightest pixels in 
the aperture, disregarding the centroid of the extraction region. 
If the emission from different wavelength bands originates from 
different spatial locations in the aperture, this will likely be done 
incorrectly for some of them (Sahnow 2010, private communi- 
cation). In Fig. A.l the danger of performing the astigmatism 
correction on the slope of the detector background in SiC IB is 
demonstrated. It is unclear whether the shape of what appears to 
be the target signal in a) is caused by a correction to either the 
target signal or the slope of the background. The SiC 2A signal 
on the other hand, does not seem to be affected in the same way. 
We need to consider that the LyC could be emitted from an off- 
center star-forming region, the uncertainty of the astigmatism 
correction, and the risk that not all the signal falls within the 
point source extraction window, hence we need to consider both 
cases (point source and extended source) in our analysis. 

A.3. From 2D to 1D 

In the analysis we also wished to investigate the reliability of 
the ID extraction routine, and tried different options such as 
optimal- or normal extraction, and selecting or not the back- 



ground subtraction routine. There are eight ID spectrum files 
produced by CalFUSE, one for each spectrum. These provide 
several options for the user to display the data, where for ex- 
ample FLUX, COUNTS, WEIGHTS, and BKGD (background) 
are offered. The WEIGHTS column represents the raw counts 
in COUNTS corrected for dead-time effects, but for faint targets 
such as Haro 11 the dead-time correction is negligible, so the 
two are essentially identical. The final flux is calculated by sub- 
tracting BKGD from WEIGHTS. 

The data were reduced with different parameter set-ups. First 
we used the extended source setting (the default used in |BergvaTI] 
et al. (2006)). This can be seen in Fig A. 2 a), where the total SiC 
IB signal (WEIGHTS) is overplotted with the fitted background 
(BKGD, dotted line). A clear overestimate of the background is 
visible at the shorter wavelengths. The same is plotted for SiC 
2A in panel e). [Dixon et al. (2007i argue for using the point 
source setting whenever possible in order to improve the back- 
ground model for faint targets. We tried this both with normal 
extraction (all pixels get equal weight) as well as with the stan- 
dard optimal extraction routine where central pixels get more 
weight (the latter being the case used in Grime s~et al.| ([2007)). 
The resulting SiC IB spectra can be seen in Fig. A. 2 b) - d) for 
IB and f) - h) for 2 A. In b) and f) the normal extraction produces 
a strong increase of the background at the shorter wavelengths. 
When the background is subtracted, this results not only in an 
unphysical negative flux, but also an artificially produced Ly- 
limit feature that makes it impossible to determine the zero-level 
of the flux below the Milky Way Ly-limit (912 A). In c) - d) and 
g) - h) the point source setting was used with the optimal ex- 
traction routine. The difference between them is that in c) and g) 
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the WEIGHTS column is taken from a file where the background 
routine was turned off. The overplotted backgrounds are instead 
taken from the files in d) and h). The WEIGHTS columns should 
be identical regardless of whether the background has been sub- 
tracted or not, thus plot c) is expected to be identical to d), and g) 
identical to h), which they clearly are not. The reason for this is 
not known but could possibly, according to Dixon (2007, private 
communication), be caused by the complex optimal extraction 
routine involving probability arrays used on low S/N spectra. 

Comparing Fig A. 2 to Fig. [8] we see a large difference be- 
tween the two background models in SiC 2A. In Fig. [8] where 
the background has been fitted directly on the detector, the curve 
is smooth. In Fig A. 2 on the other hand, the shape of the tem- 
plate background file is seen to vary in a way often uncorrelated 
to the Haro 1 1 data. This demonstrates the danger of using non- 
simultaneous template background files, and for faint targets the 
effect can be severe. The features will have a large effect on the 
output flux for which the amplitude will vary in an unknown 
way over the spectrum. The effect is largest at the shortest wave- 



lengths in all panels of Fig A. 2 
Summary Appendix A: 

We can conclude that both the point source and extended source 
options must be considered for the Haro 1 1 data since there are 
benefits and risks with both. While the point source option might 
be unable to include all of the signal, might skew the weighting 
if the emission is off-center, and in the case of SiC IB perform 
a questionable astigmatism correction, the extended source op- 
tion, which is otherwise more stable, has been shown to include 
more of the background region and therefore increase the un- 
certainty in the background model. We have also shown that 
the non-simultaneous background templates used by CalFUSE 
v3.2 display features that are uncorrelated to the Haro 11 data, 
and that these features can cause large errors in the final flux at 
random places over the spectrum (but especially at the shortest 
wavelengths). This analysis has therefore led us to develop our 
own background model from a direct fit to the detector response 
(Section |3.2| i, which is the standard procedure when reducing 
astronomical data. 



